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Owing to the enormous interest the rapidly growing field of topological states of matter (TSM) 
has attracted in recent years, the main focus of this review is on the theoretical foundations of TSM. 
Starting from the adiabatic theorem of quantum mechanics which we present from a geometrical 
perspective, the concept of TSM is introduced to distinguish gapped many body ground states of 
non-interacting systems and mean field superconductors, respectively, regarding their global geomet- 
rical features. These classifying features are topological invariants defined in terms of the adiabatic 
curvature of these bulk insulating systems. We review the general classification of TSM in all sym- 
metry classes in the framework of K-Theory. Furthermore, we outline how interactions and disorder 
can be included into the theoretical framework of TSM by reformulating the relevant topological 
invariants in terms of the single particle Green's function and by introducing twisted boundary con- 
ditions, respectively. We finally integrate the field of TSM into a broader context by distinguishing 
TSM from the concept of topological order which has been introduced to study fractional quantum 
Hall systems. 



I. INTRODUCTION 

Stimulated by the theoretical predictiorP^ and ex- 
perimental discovery^ of the quantum spin Hall (QSH) 
state, tremendous interest has been recently attracted by 
the study of topologocical properties of non-interacting 
band structures'"^. A topological state of matter (TSM) 
is a non-interacting band insulator which is topologically 
distinct from a conventional insulator. This means that 
a TSM cannot be adiabatically, i.e., while maintaining a 
finite bulk gap, deformed into a conventional insulator 
without breaking its fundamental symmetries. This 
classification approach provides a new paradigm in con- 
densed matter physics which goes beyond the mechanism 
of local order parameters associated with spontaneous 
symmetry breaking^ Different TSM differ in the value 
of their defining topological invariant but concur in all 
conventional symmetries. Interestingly, these global 
topological features are not always immediately visible 
in the microscopic equations of motion. However, the 
bulk topology leads to unique finite size effects at the 
boundary of a finite sam ple wh ich has been coined bulk 
boundary correspondence ^ 9 ^ 10 !. This general mechanism 
gives rise to peculiar holographic transport properties 
of topologically non-trivial systems. Predicting and 
probing the rich phenomenology of these topological 
boundary effects in mesoscopic samples has become one 
of the most rapidly growing fields in condensed matter 
physics in recent years. 

The focus of this review is threefold. First, we view the 
topological invariants characterizing TSM as the global 
analogues of geometric phases associated with the adi- 
abatic evolution of a physical systenfSE^. Geometric 
phases are well known to reflect the local inner-geometric 
properties of the Hilbert sp ace of a system with a param- 
eter dependent Hamiltoniar i 14 * 15 l The topological invari- 



ants represent global features of a Bloch Hamiltonian, 
i.e., of a system whose parameter space is its fc-space. 
This approach integrates TSM and physical phenomena 
related to geometric phases into a common theoretical 
context. Second, we present the topological classifica- 
tion of TSM in the framework of K-Theory^ in a more 
accessible way than in the pioneering work by Kitae\^. 
Third, we discuss very recent developments concerning 
the generalization of the relevant topological invariants 
to interacting and disordered systems several of which 
were published after existing reviews on TSM. 



II. THE ADIABATIC THEOREM 

The gist of the adiabatic theorem can be understood 
at a very intuitive level: Once prepared in an instanta- 
neous eigenstate with an eigenvalue which is separated 
from the neighboring states by a finite energy gap A, the 
system can only leave this state via a transition which 
costs a finite excitation energy A. A simple way to esti- 
mate whether such a transition is possible is to look at 
the Fourier transform "H(u;) of the time dependent Hamil- 
tonian H(t). If the time dependence of rl is made suffi- 
ciently slow, H(w) will only have finite matrix elements 
for u) <C A. In this regime the system will stick to the 
same instantaneous eigenstate. This behavior is known 
as the adiabatic assumption. 



A. Proof due to Born and Fock 

The latter rather intuitive argument is at the heart of 
the adiabatic theorem of quantum mechanics which has 
been first proven by Born and Fock in 1928^ for non- 
degenerate systems. Let {\n(t))} n be an orthonormal 
set of instantaneous eigenstates of "H(i) with eigenvalues 
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{E n (t)} n . The exact solution of the Schrodinger equation 
can be generally expressed as 



i*w> = E c «( t )K < )) c 
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where the dynamical phase 4>%(t) = f t E n (r)dT has 
been separated from the coefficients c„ (t) for later conve- 
nience. Plugging Eq. ([I]) into the Schrodinger equation 
yields 



- c n{n\j t \n) - ^ 
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with 0™" 



(*) = $d(0 — ^dW- The salient consequence 
of the adiabatic theorem is that the last term in Eq. 
([2]) can be neglected in the adiabatic limit since its de- 
nominator \E n — E m \ > A is finite whereas the ma- 
trix elements of become arbitrarily small. More 
precisely, if we represent the physical time as t = Ts, 
where s is of order 1 for a change in the Hamiltonian 
of order A and T is the large adiabatic timescale, then 
It = Tls- Now, 3j%(i(s)) is by construction of order 
A. The entire last term in Eq. |2) is thus of order ij,. 
Under these conditions^, Born and Fock^ showed that 
the contribution of this second term vanishes in the adi- 
abatic limit T — > 00. Note that this is not a trivial result 
since the differential equation (|2| is supposed to be in- 
tegrated from t — to t ~ T, so that one could naively 
expect a contribution of order 1 from a coefficient that 
scales like 1/T. The coefficient of c n in the first term on 
the right hand side of Eq. p| is purely imaginary since 

= li( n \ n ) = (s( n Dl n ) + v l IJjl n ) an d hence does not 
change the modulus of c n when the differential equation 
c„ = —Cn(n\4r\n) is solved as 



c n (t) 



Is- \ - (n\4-\n)dT 



(3) 



n\jjj\n) — Vt amounts to 



Born and FoclP^ argue that 
a choice of phase for the eigenstates and therefore ne- 
glect also the first term on the right hand side of Eq. 

This review article is mainly concerned with physical 
phenomena associated with corrections to this in general 
unjustified assumption. 



B. Notion of the geometric phase 

By the latter assumption, Ref. [18] overlooks the po- 
tentially nontrivial adiabatic evolution, known as Berry's 
phased, associated with a cyclic time dependence of H. 
After a period [0, T] of such a cyclic evolution, Eq. ^ 
yields 



c n {T) =c n (0)e-£> T W^I«>^ 



(4) 



To understand why the phase factor e" 

in general not be gauged away, we remember that 



the Hamiltonian depends on time via the time depen- 
dence R(t) of some external control parameters. Hence, 
{n\4r\n) = (n\d^\n)R^, where 9 M = gj^j. To reveal the 
mathematical structure of the latter expression, we define 



A 1 



dt 



= A^Rf i = -i(n\d M \n)R' t , 



(5) 



where A B = A B dR^ is called Berry's connection. 
A B clearly has the structure of a gauge field: Under the 
local gauge transformation \n) — > e^|n) with a smooth 
function R H> £(R), Berry's connection transforms like 

A B -> A B + dt 

Furthermore, the cyclic evolution defines a loop 7 : t 1— ¥ 
R(t), t e [0,T] , R(0) = R(T) in the parameter manifold 
1Z. 7 can be expressed as the boundary of some piece 
of surface S. Using the theorem of Stokes, we can now 
calculate 



(n\ — \n)&r 



A 1 



dA 



(6) 



where in the last step Berry's curvature T B 
dR v is defined as 



T B u dR»A 



-i {{d^n\d v n) - {d v n\d^n)) = 2lra.{(d ll n\d v n)} 



with the shorthand notation \dfjji) = d^\n). Note that 
T B is a gauge invariant quantity that is analogous to the 
field strength tensor in electrodynamics. Defining the 
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Berry phase associated with the loop 7 as ip B = J A 
J s T B we can rewrite Eq. Q as 

c n (T) = c n (0)e-<. (7) 

The manifestly gauge invariant Berry phase Lp B can have 
observable consequences due to interference effects be- 
tween coherent superpositions that undergo different adi- 
abatic evolutions. The analogue of this phenomenology 
due to an ordinary electromagnetic vector potential is 
known as the Aharonov-Bohm effectpQj. The geometri- 
cal reason why Berry's connection A B cannot be gauged 
away all the way along a cyclic adiabatic evolution is the 
same as why a vector potential cannot be gauged away 
along a closed path that encloses magnetic flux, namely 
the notion of holonomy on a curved manifold. We will 
come back to the concept of holonomy shortly from a 
more mathematical point of view. For now we only com- 
ment that the Berry phase ip B is a purely geometrical 
quantity which only depends on the inner-geometrical re- 
lation of the family of states \n (R)) along the loop 7 and 
reflects an abstract notion of curvature in Hilbert space 
which has been defined as Berry's curvature T B . 



C. Proof due to Kato 

For a degenerate eigenvalue, Berry's phase is promoted 
to a unitary matrix acting on the corresponding degener- 
ate eigenspace^l. The first proof of the adiabatic theorem 
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of quantum mechanics that overcomes both the limita- 
tion to non-degenerate Hamiltonians and the assumption 
of an explicit phase gauge for the instantaneous eigcn- 
states was reported in the seminal work by Tosio KatcW 
in 1950. We will review Kato's results briefly for the 
reader's convenience and use his ideas to illustrate the 
geometrical origin of the adiabatic phase. The explicit 
proofs are presented at a very elementary and self con- 
tained level in Ref. 11. Our notation follows Ref. IHH 
which is convenient to relate the physical quantities to 
elementary concepts of differential geometry. 

Let us assume without loss of generality that the sys- 
tem is at time to = in its instantaneous ground state 
|\&o(0)) or, more generally, since the ground state might 
be degenerate, in a state |$) satisfying 



P(0)|*) = |*), 



(8) 



where P(t) is the projector onto the eigenspace associated 
with the instantaneous ground state energy Eo(t) which 
is defined as 



P(t) 



dz 



1 



2m J c z-H(t) 



where the complex contour c encloses Eq(£) which is again 
assumed to be separated from the spectrum of excita- 
tions by a finite energy gap A > 0. To understand the 
adiabatic evolution, we are not interested in the dynam- 
ical phase 4>D{t) — JoE>o(T)dT. We thus define a new 
time evolution operator U(t,0) = e i<t>D( ^U{%0) . Clearly, 
Li represents the exact time evolution operator of a sys- 
tem which has the same eigenstates as the original system 
but has been subjected to a time dependent energy shift 
that transforms Eo(t) — > Eo(t) = V t . Kato proved the 
adiabatic theorem in a very constructive way by writing 
down explicitly the generator A of the adiabatic evolu- 
tion: 



P,P 



(9) 



In the adiabatic limit, U(t,0)P(0) was shown W to con- 
verge against the adiabatic Kato propagator /C, i.e., 



W(t,0)P(0) 



adiabatic 



V^K^OHTe-^-^K (10) 



The adia batic assumption is now a direct corollary from 
Eq. ( 10 ) and can be elegantly expressed aJ^H 



P(t)JC(t,0) =£(i,0)P(0), 



(11) 



implying that a system, which is prepared in an instan- 
taneous ground state at to — 0, will be propagated to a 
state in the subspace of instantaneous ground states at 
t by virtue of Kato's propagator K. Note that K, is a 
completely gauge invariant quantity, i.e., independent of 
the choice of basis in the possibly degenerate subspace of 
ground states. The Kato propagator JC(T, 0) associated 
with a cyclic evolution in parameter space thus yields the 



Berry phased and its non-Abelian generalizatiorl^, re- 
spectively. We will call this general adiabatic phase the 
geometric phase (GP) in the following. The GP /C 7 rep- 
resenting the adiabatic evolution along a loop 7 in pa- 
rameter space can be expressed in a manifestly gauge 
invariant way as 



/C 7 = Te 



(12) 

Kato's propagator is the solution of an adiabatic ana- 
logue of the Schrodinger equation, an adiabatic equation 
of motion that can be written as 



h A I — 

dt \dt 



!*(*)} = 0, 



(13) 



for states satisfying P(t)\^(t)) = \^>{t)) : i.e., states 
in the subspace of instantaneous groundstates. Before 
closing the section, we give a general and at least numer- 
ically always viable recipe to calculate the Kato propaga- 
tor K.(t, 0). We first discretize the time interval [0, t] into 



n steps by defining tj 



The discrete version of Eq. 



(131 for the Kato propagator reads (see Eq. (19 



K{U,0) - £(ii_i,0) - {(P(*i) - P(ti-i)) P(ii-i)- 
P(ti) (P(ti) - P{U_ x ) )} £(«,_ x , 0) . (14) 

Using P(ti-i)JC(ti-i,Q) = /C(tf_i,0) and P 2 = P, Eq. 
( 14 1 can be simplified to 



K(ti,0) = P(t i )/C(t i _i,0), 

which is readily solved by K (ij, 0) = Y[ l j =0 P(tj)- Taking 
the continuum limit yield j^- 14 * 21 ! 



K(t,0) = lim TTPfe), 

77, — J-OO A 



(15) 



i=0 



which is a valuable formula for the practical calculation 
of the Kato propagator. 



III. GEOMETRIC INTERPRETATION OF 
ADIABATIC PHASES 

We now view the adiabatic time evolution as an ab- 
stract notion of parallel transport in Hilbert space and re- 
veal the GP associated with a cyclic evolution as the phe- 
nomenon of holonomy due to the presence of curvature 
in the vector bundle of ground state subspaces over the 
manifold 1Z of control parameters. Interestingly, Kato's 
approach to the problem provides a gauge invariant, i.e., 
a global definition of the geometrical entities connection 
and curvature, whereas standard gauge theories are de- 
fined in terms of a complete set of local gauge fields along 
with their transition functions defined in the overlap of 
their domains. This difference has an interesting physi- 
cal ramification: Quantities that are gauge dependent in 
an ordinary gauge theory like quantum chromodynamics 



4 



(QCD) are physical observables in the theory of adia- 
batic time evolution. To name a concrete example, only 
gauge invariant quantities like the trace of the holonomy, 
also known as the Wilson loop, are observable in QCD 
whereas the holonomy itself, in other words the GP de- 
fined in Eq.( 12 ), is a physical observable in Kato's theory. 



This subtle difference h as be en overlooked in standard 
literature on this subjedP^H which we interpreted as an 
incentive to clarify this point below in greater detail. 



A. Adiabatic time evolution and parallel transport 

To get accustomed to parallel transport, we first 
explain the general concept with the help of a very 
elementary example, namely a smooth piece of two 
dimensional surface embedded in M 3 . If the surface is 
flat, there is a trivial notion of parallel transport of 
tangent vectors, namely shifting the same vector in the 
embedding space from one point to another. However, 
on a curved surface, this program is ill-defined, since 
a tangent vector at one point might be the normal 
vector at another point of the surface. Put shortly, 
a tangent vector can only be transported as parallel 
as the curvature of the surface admits. On a curved 
surface, parallel transport along a curve is thus defined 
as a vanishing in-plane component of the directional 
derivative, i.e., a vanishing covariant derivative of a 
vector field along a curve. The normal component of 
the directional derivative reflects the rotation of the 
entire tangent plane in the embedding space and is 
not an inner-geometric quantity of the surface as a two 
dimensional manifold. 



The analogue of the curved surface in the context of 
adiabatic time evolution is the manifold of control pa- 
rameters 1Z, parameterizing for example external mag- 
netic and electric fields. The analogue of the tangent 
plane at each point of the surface is the subspace of de- 
generate ground states of the Hamiltonian H(R) at each 
point R in parameter space. An adiabatic time depen- 
dence of H. amounts to traversing a curve t \— > R(t) in 
1Z at adiabatically slow velocity. A cyclic evolution is 
uniquely associated with a loop 7 in 1Z. We will now ex- 
plicitly show that the adiabatic equation of motion ( 13 ) 



defines a notion of parallel transport in the fiber bundle of 
ground state subspaces over 72. in a completely analogous 
way as the ordinary covariant derivative V on a smooth 
surface defines parallel transport in the tangent bundle of 
the smooth surface. We first note that 4? — R^d^ is re- 
ferring to a particular direction i? M in parameter space, 
which depends on the choice of the adiabatic time de- 
pendence of %. We can get rid of this dependence by 
rephrasing Eq. ( 13 ) as 



(d + A) |*> = 0, 



(16) 



for notational convenience. The adiabatic derivative 
D = d + A takes a tangent vector, e.g., J^, as an argu- 
ment to boil down to the directional adiabatic derivative 
A 1 A (A 

dt ' V dt 



appearing in Eq. (13 1. For the following 
analysis the identities P 2 = P and P\^S) — |^) are of key 
importance. It is now elementary algebra to show 



P{dP)P = 0. 



(17) 



Eq. (17) has a simple analogue in elementary geometry: 



Consider the family of unit vectors {n(t)} t where t pa- 
rameterizes a curve on a smooth surface. Then, since 
1 = (n|n), we get = j^{n\n) = 2{n\h), i.e., the change 
of a unit vector is perpendicular to the unit vector itself. 
Using Eq. (17), we immediately derive PAl'i') = and 
with that 



D\V) = & Pd\^) = 0. 



(18) 



This makes the analogy of our adiabatic derivative D = 
d + A to the ordinary notion of parallel transport mani- 
fest: 1^) is parallel-transported if the in-plane component 
of its derivative vanishes. 



Curvature and holonomy 

Let us again start with a very simple example of a 
curved manifold, a two dimensional sphere S 2 , which 
has constant Gaussian curvature. Parallel-transporting 
a tangent vector around a geodesic triangle, say the 
boundary of an octant of the sphere gives a defect angle 
which is proportional to the area of the triangle or, more 
precisely, the integral of the Gaussian curvature over the 
enclosed area. This defect angle is called the holonomy 
of the traversed closed path. This elementary example 
suggests that the presence of curvature is in some sense 
probed by the concept of holonomy. This intuition is 
absolutely right. As a matter of fact, the generalized 
curvature at a given point x of the base manifold of 
a fiber bundle is defined as the holonomy associated 
with an infinitesimal loop at x. More concretely, the 
curvature fl is usually defined as Q,^ v = [V^V^] which 
represents an infinitesimal parallel transport around a 
parallelogram in the /i^-plane. 

In total analogy, we define 

•7^1*) - [D ll ,D v ]\*)=P[P ll ,P v ]P\*), (19) 

with the shorthand notation P^ = d^P. Restricting the 
domain of T to states which are in the projection P, we 
can rewrite Eq. ( 19 1 as the operator identity 



where A = —[(dP),P] and here as in the following 
P|\]/) = J if) and the ^-dependence has been dropped 



T = T^dR^ A dR v — P [(dP), {dP)} P, (20) 

where the product of the two differential forms dP in the 
commutator is to be understood as the usual exterior 
A-product. 
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In the general case of a non-Abelian adiabatic connec- 
tion, i.e., if the dimension of P is larger than 1, we cannot 
simply use Stokes theorem to reduce the evaluation of Eq. 
(12 1 to a surface integral of T over the surface bounded 



Eq. (24 1 makes the relation between Kato's formulation 



by 7, as has been done in the case of the Abelian Berry 
curvature in Eq. ([6]). However, the global one to one cor- 
respondence between curvature and holonomy still exists 
and is the subject of the Ambrose-Singer theorem^. 



Relation between Kato's and Berry's language 

In order to make contact to the more standard lan- 
guage of gauge theory, we will now express Kato's mani- 
festly gauge invariant formulation^ in loc al coo rdinates 
thereby recovering Berry's connection ,A^ 12 I 14 I anc l its 
non-Abelian generalization^, respectively. For this pur- 
pose, let us fix a concrete basis {\a(R))} a , R € O C 1Z in 
an open subset O of the parameter manifold. We assume 
the loop 7 to lie inside of O. Otherwise we would have 
to switch the gauge while traversing the loop. We will 
drop the independence of \a) right away for notational 
convenience. The projector P can then be represented as 
P = J2 a \ a )( a \- Let us start the cyclic evolution without 
loss of generality with |W(0)) = |cc(0)). From Eq. (11) we 
know that the solution = JC(t, 0)|a(0)) of Eq. Q 

satisfies P(t)\^>(t)) = \$f(t)) at every point in time during 
the cyclic evolution. Hence, we can represent ^(t)) in 
our gauge as 



|*>=£O3|*>I/J> = 0) 





(21) 



where the i^ependence has been dropped for brevity. 
From Eq. (18 1, we know that PAXQ) — which implies 



(7|4|*) = 0. Plugging this into Eq. (21) yields 



±UB 
dt la 



£<7l 



(22) 



Redefining A B for the non-Abelian case as a matrix val- 
ued gauge field through A B p = -i (a \d^\P)dR^, Eq. Q 
is readily solved as 

U B (t)=Te-*fo AB ^ dT . 

The representation matrix of the GP associated with the 
loop 7 then reads 



U» = Te 



(23) 



By construction, U B is the representation matrix of the 
GP /C 7 , i.e., 

(U*) ap = (a(0)|/C 7 |/3(0)), 
or, more general, for any point in time along the path 
(U B (t)) a0 = {a(t)\lC{t, O)|0(O)>. (24) 



of adiabatic time evolution and the non-Abelian Berry 
phase manifest. In contrast to the gauge independence 
of Kato's global connection A, A B behaves like a local 
connection (see Ref. [M] for rigorous mathematical def- 
initions) and depends on the gauge, i.e., on our choice 
of the family {\a{R))} a of basis states. Under a smooth 
family of basis transformations {£^(i?Hg acting on the 
local coordinates A B transforms lik e! 23 ! 24 ! 

A B ^ A B = U~ 1 A B U + U~ 1 dU (25) 

resulting in the following gauge dependence of Eq. (p3| , 



U B -+ U B 



U~ l U B U, 



(26) 



which only depends on the basis choice U = U (R(0)) at 
the starting point of the loop 7. 

Inserting our representation P = Yla\ a )( a \ into the 
gauge independent form of the curvature, Eq.(20), we 
readily derive 



T, 



B,a(3 _ 



„., _ x a\ [P„ P v ] \P) = (dA B )f„ + (A B A A B )f v , 



which defines T B as the usual curvature of a non-Abelian 
gauge fielcP^l, i.e., 



F B = dA B + A A A 



(27) 



which transforms under a local gauge transformation 
U like 

J~ B -> XJ- X T B XJ. 



B. Gauge dependence and physical observability 

The gauge depe nden ce of the non-Abelia n Ber ry 
phase U B (see Eq. (26)) has led several authors^SIlIl to 
the conclusion that only gauge independent features like 
the trace and the determinant of U B can have physical 
meaning. However, working with Kato's manifestly 
gauge invariant formulation, it is understood that the 
entire GP /C 7 is experimentally observable. In the 



remainder of this section we will try to shed some light 
on this ostensible controversy. 

In gauge theory, it goes without saying that explicitly 
gauge dependent phenomena are not immediately physi- 
cally observable and that only the gauge invariant infor- 
mation resulting from a calculation performed in a spe- 
cial gauge can be of physical significance. At a formal 
level this is a direct consequence of the fact that the La- 
grangian of a gauge theory is constructed in a manifestly 
gauge invariant way by tracing over the gauge space in- 
dices. The physical reason for this is quite simple: A 
concrete gauge amounts to a local choice of the coordi- 
nate system in the gauge space. Under a local change 
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of basis, a non-abelian gauge field A transforms like (see 
also Eq. p5}) 



A -> A = LT^AE/ + f7 -1 df7, 

where [7 (x) is a smooth family of basis transformations, 
with x labeling points in the base space of the theory, 
e.g., in Minkowski space. Now, since the gauge space is 
an internal degree of freedom, the basis vectors in this 
space are not associated with physical observables. This 
situation is fundamentally changed in Kato's adiabatic 
analogue of a gauge theory. Here, the non-Abelian struc- 
ture is associated with a degeneracy of the Hamiltonian, 
e.g., Kramers degeneracy in the presence of time rever- 
sal symmetry (TRS). For a system in which spin is a 
good quantum number, Kramers degeneracy is just spin 
degeneracy, which makes the spin the analogue of the 
gauge degree of freedom in an ordinary gauge theory. 
However, the magnetic moment associated with a spin is 
a physical observable which can be measured. The basis 
vectors, e.g., |t), ||) have an objective meaning for the ex- 
perimentalist (a magnetic moment that points from the 
lab-floor to the sky which we call z- direction). For con- 
creteness, let us assume that we have calculated a GP 
^-7 = It) (II + 11)01- The representation matrix of /C 7 in 
this basis of S z eigenstates is clearly the Pauli matrix 
a x . Choosing a different gauge, i.e., a different basis for 
the gauge degree of freedom at the starting point of the 
cyclic adiabatic evolution, we of course would have ob- 
tained a different representation matrix for K.^, e.g., 
cr 2 , had we chosen the basis as eigenstates of S x (see Eq. 



(26)). However, the fact that K. 1 rotates a spin which is 



initially pointing to the lab-ceiling upside down is gauge 
independent physical reality. 



IV. FROM GEOMETRY TO TOPOLOGY 

In Section |III| we worked out the relation between 
the GP and the notion of curvature as a local geometric 
quantity. The topological invariants characterizing TSM 
are in some sense global GPs. They measure global 
properties which cannot be altered by virtue of local 
continuous changes of the physical system. Continuous is 
at this stage of the analysis synonymous with adiabatic, 
i.e., happening at energies below the bulk gap. Later 
on, we will additionally require local continuous changes 
to respect the fundamental symmetries of the physical 
system, e.g., particle hole symmetry (PHS) or time 
reversal symmetry (TRS). 

Let us illustrate the correspondence between local cur- 
vature and global topology of a manifold with the help 
of the simplest possible example. We consider a two di- 
mensional sphere S 2 with radius r. This manifold has a 
constant Gaussian curvature of n = The integral of 
k over the entire sphere obviously gives 47T, independent 
of r. The Gauss-Bonnet theorem in its classical form 
(see, e.g., Ref. |2^|) relates precisely this integral of the 



Gaussian curvature of a closed smooth two dimensional 
manifold M. to its Euler characteristic \ m the following 
way: 



1 

27 



M 



X(M) 



(28) 



Note that \ is a purely algebraic quantity which is de- 
fined as the number of vertices minus the number of edges 
plus the number of faces of a triangulatiorPS of M. . x is 
by construction of simplicial homologjEH a topological in- 
variant which can only be changed by poking holes into 
M. and gluing the resulting boundaries together so as 
to create closed manifolds with different genus. Hence, 
Eq. ( 28 1 nicely demonstrates how the integral of the lo- 



cal inner-geometric quantity k over the entire manifold 
yields a topological invariant. Concretely, for our exam- 
ple S 2 , a triangulation is provided by continuously de- 
forming the sphere into a tetrahedron. Simple counting 
of vertices, edges, and faces yields x(S 2 ) = 4 — 6 + 4 = 2, 
in agreement with Eq. ( 28 1 . More generally speaking, 
e = £- is our first encounter with a characteristic clasiP-H 
the so called Euler class of Ai, which upon integration 
over M. yields the topological invariant \- Similar math- 
ematical structures will be ubiquitous when it comes to 
the classification of TSM. The simplest example in this 
context is the quantum anomalous Hall (QAH) stated, 
a 2D insulating state in symmetry class A (see Section 
VA) which is characterized by its first Chern number 



see Section VI A for a detailed discussion) 



BZ 



iTr [JF] 
2?r 



(29) 



The formal analogy between Eq. (28 



and Eq. (29) is 



lich the adiabatic 



(29) is the fc-space of 



striking. The parameter space over w 
curvature T is integrated in Eq. 
the physical system, i.e., the Brillouin zone (BZ) for a 
periodic system, which has the topology of a torus. In 
Section [IlJ we showed that the GP can be viewed as the 
flux through the surface in parameter space which is 
bounded by the corresponding cyclic adiabatic evolution. 
Along similar lines, the first Chern number is analogous 
to a monopole charge enclosed by the entire BZ of the 
QAH insulator. 

The close correspondence between adiabatic evolution 
and the first Chern number has also been discussed in 
RefspjHUl p or a 2j) S y S tem in cylinder geometry, a cir- 
cumferential electric field can be modeled by a time de- 
pendent magnetic flux threading the cylinder in axial di- 
rection. In this scenario, the first Chern number can 
be viewed as the shift of the charge polarization in ax- 
ial direction associated with the adiabatic threading of 
one flux quantum. Along these lines, LaughliiP^ had in- 
terpreted the quantum Hall effect as an adiabatic charge 
pumping process, even before the formal relation between 
the Hall conductivity and the first Chern number was es- 
tablished in RefsPtSl. 
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V. BULK CLASSIFICATION OF ALL 
NON-INTERACTING TSM 

Very generally speaking, the understanding of TSM 
can be divided into two subproblems. First, finding 
the group that represents the topological invariant for 
a class of systems characterized by their fundamental 
symmetries and spatial dimension. Second, assigning 
the value of the topological invariant to a representative 
of such a symmetry class, i.e., measuring to which 
topological equivalence class a given system belongs. We 
address the first problem in this section and the second 
problem in |VI| 



The general idea that yields the entire table of TSM 
is quite simple: In addition to requiring a bulk insulat- 
ing gap, physical systems of a given spatial dimension 
are divided into 10 symmetry classes distinguished by 
their fundamental symmetries, i.e., TRS, PHS, and chiral 
symmetry (CSJPH The topological properties of the cor- 
responding Cartan symmetric spaces of quadratic candi- 
date Hamiltonians determine the group of possible topo- 
logically inequivalent systems. We outline the mathemat- 
ical structure behind this general classification scheme in 
some detail. First, we briefly review the construction of 
the ten universality classes^. Then, we present the asso- 
ciated topological invariants for non-interacting systems 
of arbitrary spatial dimension giving a complete list of all 
TSM 6 that can be distinguished by virtue of this frame- 
work. Finally, we discuss in some detail the origin of 
characteristic patterns appearing in this table using the 
framework of K-Theory along the lines of the pioneering 
work by Kitaev^. 



A. Cartan-Altland-Zirnbauer symmetry classes 

A physical system can have different types of symme- 
tries. An ordinary symmetry-Sis characterized by a set of 
unitary operators representing the symmetry operations 
that commute with the Hamiltonian. The influence of 
such a symmetry on the topological classification can 
be eliminated by transforming the Hamiltonian into a 
block-diagonal form with symmetry-less blocks. The 
total system then consists of several uncoupled copies 
of symmetry-irreducible subsystems which can be clas- 
sified individually. In contrast, the "extremely generic 
symmetries'^ follow from the anti-unitary operations of 
TRS and PHS. Involving complex conjugation according 
to Wigner, they impose certain reality conditions on 
the system Hamiltonian. In total, the behavior of the 
system under these operations, and their combination, 
the CS operation, defines ten universality classes which 
we call the Cartan-Altland-Zirnbauer (CAZ) classes. 
For disordered systems, these classes correspond to 
ten distinct rcnormalization group (RG) low energy 
fixed points in random matrix theory^. The spaces of 
candidate Hamiltonians within these symmetry classes 



correspond to the ten symmetric spaces introduced by 
Cartan in 1926 36 defined in terms of quotients of Lie 
groups represented in the Hilbert space of the system. 
For translation-invariant systems, the imposed reality 
conditions are inherited by the Bloch Hamiltonian h(k). 

In the following, the anti-unitary TRS operation will 
be denoted by T and the anti-unitary PHS will be de- 
noted by C. The Hamiltonian H of a physical system 
satisfies these symmetries if 



THT^ 1 = H, 



and 



CHC- 



-n, 



(30) 



(31) 



respectively. According to Wigner's theorem, these anti- 
unitary symmetries can be represented as a unitary op- 
eration times the complex conjugation K. We define 
T = TK, C = CK. Using the umtarity of T, C along 
with T-L = Ti) we can rephrase Eqs. ( 30pl ) as 



TH T T* = H, 

cn T c^ = -n. 



(32) 



There are two inequivalent realizations of these anti- 
unitary operations distinguished by their square which 
can be plus identity or minus identity. For example, 
T 2 = ±1 for the unfolding of a particle with integer/half- 
integer spin, respectively. Clearly, T 2 = ±1 TT* 
■ 1 and C 2 = ±1 CC* = ±1. In total, there are thus 
nine possible ways for a system to behave under the two 
anti-unitary symmetries: each symmetry can be absent, 
or present with square plus or minus identity. For eight 
of these nine combinations, the behavior under the com- 
bination TC is fixed. The only exception is the so called 
unitary class which breaks both PHS and TRS and can 
either obey or break their combination, the CS. This class 
hence splits into two universality classes which add up to 
a grand total of ten classes shown in Tab. [TJ For a pe- 
riodic system, symmetry constraints similar to Eq. (32) 
hold for the Bloch Hamiltonian h(k), namely 



Th T {-k)T ] = h(k), 
Ch T (-k)C^ = -h(k), 



(33) 



where T, C now denote the representation of the unitary 
part of the anti-unitary operations in band space. 

For a continuum model, the real space Hamiltonian 
H{x) is defined through 



H 



d d x&(x)H(x)$>(x), 



where \& is a vector/spinor comprising all internal degrees 
like spin, particle species, etc.. The fc-space on which the 
Fourier transform H(k) of H(x) is defined does not have 



Class 


TRS 


PHS 


cs 


A (Unitary) 











AI (Orthogonal) 


+1 








All (Symplectic) 


-1 








AIII (Chiral Unitary) 








1 


BDI (Chiral Orthogonal) 


+1 


+1 


1 


CII (Chiral Symplectic) 


-1 


-1 


1 


D 





+1 





C 





-1 





Dili 


-1 


+1 


1 


CI 


+1 


-1 


1 



TABLE I. Table of the CAZ universality classes. denotes 
the absence of a symmetry. For PHS and TRS, ±1 denotes the 
square of a present symmetry, the presence of CS is denoted 
by 1. The last four classes are Bogoliubov deGennes classes 
of mean field superconductors where the superconducting gap 
plays the role of the insulating gap. 



the topology of a torus like the BZ of a periodic system. 
However, the continuum models one is concerned with in 
condensed matter physics are effective low energy/large 
distance theories. For large k, H(k) will thus generically 
have a trivial structure so that the fc-space can be en- 
dowed with the topology of the sphere S d by a one point 
compactification which maps k — > oo to a single point. 
The symmetry constraints on H(k) have the same form 
as those on the Bloch Hamiltonian h{k) shown in Eq. 
( 33 ). By abuse of notation, we will denote both H(k) and 
hjk) by h(k). Nevertheless, we will point out several dif- 
ferences between periodic systems and continuum models 
along the way. 



B. Definition of the classification problem for 
continuum models and periodic systems 

For translation-invariant insulating systems with n oc- 
cupied and m empty bands and continuum models with 
n occupied and m empty fermion species, respectively, 
the projection P(k) = J2a=i\ u a(k)) (u a (k)\ onto the oc- 
cupied states is the relevant quantity for the topological 
classification. The spectrum of the system is not of in- 
terest for adiabatic quantities as long as a bulk gap be- 
tween the empty and the occupied states is maintained. 
We thus deform the system adiabatically into a flat band 
insulator, i.e., a system with eigenenergy e_ = — 1 for all 
occupied states and eigenenergy e + = +1 for all empty 
states. The eigenstates are not changed during this defor- 
mation. The Hamiltonian of this flat band system then 
reads™! 

Q(k) = (+1) (1 - P(k)) + (-l)P(k) = 1 - 2P(k) 

Obviously, Q 2 — 1, Tr [Q] = m — n. Without further 
symmetry constraints, Q is an arbitrary U(n + m) matrix 
which is defined up to a U(n) x U(m) gauge degree of 



freedom corresponding to basis transformations within 
the subspaces of empty and occupied states, respectively. 
Thus, Q is in the symmetric space 



G 



n+m,m 



(C) = G 



n+m,n \ 



U(n + m)/(U(n) x U(m)). 



Geometrically, the complex Grassmannian Gk,i{C) is a 
generalization of the complex projective plane and is de- 
fined as the set of /-dimensional planes through the ori- 
gin of C k . The set of topologically different translation- 
invariant insulators is then given by the group g of homo- 
topically inequivalent maps k H> Q(k) from the BZ T d of 
a system of spatial dimension d to the space G n + mjm (C) 
of possible Bloch Hamiltonians. For continuum models 
T d is replaced by S d and g is by definition given by 



ir d (G 



n-\-m,m 



(C)), 



(34) 



where the n-th homotopy group ir n of a space is by 
definition the group of homotopically inequivalent 
maps from S d to this space. For translation-invariant 
systems defined on a BZ, the classification can be 
more complicated than Eq. ( 34 1 if the lower homo- 
topy groups ir s , s = 1, ...,d— 1 are nontrivial. For 
the quantum anomalous Hall (QAH) insulator^, a 
2D translation-invariant state which does not obey 
any fundamental symmetries, we can infer from 
7T2 (Gn+m,m(&)) = (G„ +m , m (C)) = that an inte- 

ger topological invariant must distinguish possible states 
of matter in this symmetry class, i.e., possible maps 
T 2 G n+miTn (C). The condition n 1 (G n+m , im (C)) = is 
necessary because the 7T2 classifies maps from S 2 , which 
is only equivalent to the classification of physical maps 
from T 2 if the fundamental group tt\ of the target space 
is trivial. The difference between the base space of a 
periodic systems which is a torus and of continuum 
models which has the topology of a sphere has interesting 
physical ramifications: The so called weak topological 
insulators are only topologically distinct over a torus 
but not over a sphere. Physically, this is visible in the 
lacking robustness of these TSM which break down with 
the breaking of translation symmetry. 

Requiring further symmetries as appropriate for the 
other nine CAZ universality classes is tantamount to 
imposing symmetry constraints on the allowed maps 
T d — > G n + TOim (C), k i— > Q(k) for translation- invariant 
systems and S d — > G„+ m , m (C), k i-> Q{k) for continuum 
models, respectively. The set of topologically distinct 
physical systems is then still given by the set of homo- 
topically inequivalent maps within this restricted space, 
i.e., the space of maps which cannot be continuously de- 
formed into each other without breaking a symmetry con- 
straint. For example, for the chiral classes characterized 
by CS = 1, Q can be brought into the off diagonal forrrPSl 



Q 
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with qq^ = 1, which reduces the corresponding tar- 
get space to U(n). For the chiral unitary class AIII 
without further symmetry constraints, the calculation 
of g amounts to calculating g — ltd (U(n)) — Z for 
odd d and g = na (U(n)) — for even d, respectively. 
Additional symmetries will again impose additional 
constraints on the map T d — > U(n), k q(k). 

This procedure rigorously defines the group g of topo- 
logical equivalence classes for non-interacting translation- 
invariant insulators in arbitrary spatial dimension and 
CAZ universality class. However, the practical calcula- 
tion of g can be highly non-trivial and has been achieved 
for continuum models via various subtle detours, for ex- 
ample the investigation of surface nonlinear er-models, 
in RefsP^. In the following, we outline a mathemati- 
cal brute force solution to the classification problem in 
terms of K-Theory which has been ori gina lly introduced 
in the seminal work by Kitaev in 200S 17 . This method 
will naturally explain the emergence of weak topological 
insulators. We will not assume any prior knowledge on 
K-Theory. In Tab. [TTJ we summarize the resulting groups 
of topological sectors g for all possible systems^. We no- 
tice an interesting diagonal pattern relating subsequent 
symmetry classes to neighboring spatial dimensions. Fur- 
thermore, the pattern of the two unitary classes shows a 
periodicity of two in the spatial dimension. If we had 
shown the classification for higher spatial dimensions, we 
would have observed a periodicity of eight in the spatial 
dimension for the eight real classes. The following dis- 
cussion is dedicated to provide a deeper understanding 
of these fundamental patterns as pioneered in Ref. 1171 
The m ention ed periodicities have first been pointed out 
in Refej™. 



Topological classification of unitary vector 
bundles 



In order to prepare the reader for the application of 
K-Theory and motivate its usefulness, we first formulate 
the classification problem in the language of fiber bun- 
dles. The mathematical structure of a non-interacting 
insulator of spatial dimension d is that of a vector bundle 
E M . The base manifold M of E is the d dimensional 
/c-space of the system, and the fiber over a point k given 
by the (projective) space of occupied states P(k). The 
gauge group of the bundle is U (n) , where n is the dimen- 
sion of P. If no further symmetry conditions are imposed 
(see class A in Tab. [I]), the question of how many topo- 
logically distinct insulators in a given dimension exist is 
tantamount to asking how many homotopically different 
U(n) vector bundles can be constructed over A4. This 
question can be formally answered for arbitrary smooth 
manifolds M. as we will outline now. The general idea 

is the following. There is a universal bundle £ X into 
which every bundle E can be embedded through a bundle 



ifP f : E ( such that 



(35) 



where / : Ai — > X is the map between the base manifolds 
associated with the bundle map /. That is to say every 
bundle can be represented as a pullback bundle /*^Hof 
the universal bundle £ by virtue of a suitable bundle map 
/. The key point is now that homotopically different bun- 
dles E are distinguished by homotopically distinct maps 
/. Thus, the set of different TSM is the set n [M, X] of 
homotopy classes of maps from the fc-space M. to the base 
manifold X of the universal bundle. X is also called the 
classifying space of U(n) and is given by the Grassmanian 
G N , n (C) = U(N)/(U(n) x U(N-n)) for sufficiently large 
N, i.e., N > rf + n]. To be generic in the dimension of 
the system d, we take the inductive limit X — G n (C°° ) — 
lim^oo G N , n (C) = limAr^ U(N)/(U(n) x U(N - »)). 
We thus found for the set Vect n (M,C) of inequivalent 
U(n) bundles over M. the expression 

Vect n (M,C) = ir[M,G n (C°°)}, 

which is known for some rather simple base manifolds. 
In particular for spheres S , there is a trick to calcu- 
late Vect n (S d ,C): S d can always be decomposed into 
two hemispheres which are individually trivial. The ho- 
motopy of a bundle over S d is thus determined by the 
clutching function f c defined in the overlap S^ -1 of the 
two hemispheres, i.e., along the equator of S d . Phys- 
ically, f c translates a local gauge choice on the upper 
hemisphere into a local gauge choice on the lower hemi- 
sphere and is thus a function f c : S^ 1 — > U(n). The 
group of homotopy classes of such functions is by defini- 
tion 7Td-i(U(n)). Interestingly, for n > these groups 



are given by 
Vect„(5 d ,C) = ^_ 1 (C/(n)) = 



Z, d — 1 odd 
{0} , d — 1 even 



(36) 



This periodicity of two in {d — 1) is known as the 
complex Bott periodicity. The physical meaning of Eq. 
(361 is the following: In the unitary universality class A, 



there is an integer topological invariant in even spatial 
dimension (e.g., QAH in d = 2) and no TSM in odd 
spatial dimension. 

This classification has two shortcomings. First, it can- 
not be readily generalized to other CAZ classes at this 
simple level. Second, only systems with the same number 
of occupied bands n can be compared. However, adding 
some trivial bands to the system should yield a system 
in the same equivalence class. Both sho rtcom ings can be 
overcome m the framework of K-TheorjPE2l. 



D. K-Theory approach to a complete classification 



K-Theor ^ 16 l 40 l is concerned with vector bundles which 
have a "sufficiently large" fiber dimension. This means, 



10 



Class 


constraint 


d=l 


d= 2 


d = 3 


d = 4 


A 


none 





z 





z 


AIII 


none on q 


z 





z 





AI 


Q T (k) = Q(-k) 











z 


BDI 


q*(k) = q(-k) 


z 











D 


T x Q T (k)r x = —Q(—k), m = n 


z 2 


z 








Dili 


q(k) T = —q(—k), m = n even 


z 2 


z 2 


z 





All 


ia y Q T ' (k)(—i<j v ) = Q(—k), m, n even 





z 2 


z 2 


z 


CII 


ia y q* (k)(—ia y ) = q(—k), m = n even 


z 





z 2 


z 2 


C 


T y Q T (k)r y = —Q(—k), m = n 





z 





z 2 


CI 


q(k) T — q(—k), m = n 








z 






TABLE II. Table of all groups g of topological equivalence classes. The first column denotes the CAZ symmetry class, divided 
into two unitary classes without anti-unitary symmetry (top) and eight "real" classes with at least one anti-unitary symmetry 
(bottom). The second column shows the symmetry constraints on the fiat band maps, where we have chosen the representation 
T = K for T 2 = 1, T = io y K for T 2 = — 1, as well as C — r x K for C 2 = 1, C — r y K for C 2 = —1. Here, a y denotes the 
Pauli matrix in spin space, t x , r y denote Pauli matrices in the particle hole pseudo spin space of Bogoliubov deGennes Hilbert 
spaces. In the last four columns, g is listed for d = 1, . . . , 4. 



that topological defects which can be unwound by just 
increasing the fiber dimension are not visible in the 
resulting classification scheme. This is physically reason- 
able, as trivial occupied bands from inner localized shells 
for example increase the number of bands as compared 
to the effective low energy models under investigation. 
Models of different number of such trivial bands should 
be comparable in a robust classification scheme. The 
use of K- Theory for the classification of TSM has been 
pioneered in RefsP^l and more systematically been 
discussed in Ref. [TTj 



This is because for vector bundles on a smooth manifold 
every bundle can be augmented to a trivial bundle, i.e., 



V G 3 H> , G® H = Z 



(39) 



Eq. ( 38 1 naturally leads to the notion of a subtraction 
on VkJTA) by virtue of the Grothendieck construction: 
Consider the pairs (E 1 ,E 2 ) € Vk(M) x Vk(M) and de- 
fine the equivalence relation 

(E U E 2 ) ~ (F 1 ,F 2 ) & 3 H F X ®E 2 ®H ~E 1 @F 2 @H. 

(40) 



Crash-course in K-Theory 

The direct sum of two vector bundles E(BF is the direct 
sum of their fibers over each point. This addition has only 
a semi-group structure, since E@G = F®G=f>E~F. 
A minimal counterexample is given by E = TS 2 , F = 



S' 2 



F is clearly trivial, whereas E, the tangent 



bundle of S is well known to be non-trivial. However, 
adding NS 2 , the bundle of normal vectors to S 2 , to both 
bundles E, F we obtain the same trivial bundle S 2 xl 3 . 
This motivates the concept of stable equivalence 



E~F^E< 



Ft 



(37) 



where Z n = M x K n , K = M, C is the trivial bundle over 
the fixed base manifold A4 , which plays the role of an ad- 
ditive zero as far as stable equivalence is concerned. We 
denote the set of A'-vector bundles over M by Vk(M) in 
the following. Note that stably equivalent bundles can 
have different fiber dimension, as m ^ n in general in 
Eq. (37). The benefit of this construction is: 



E®G=F®G^E~F 



(38) 



Looking at Eq. (|40j), we can intuitively think of the 
equivalence class (E 1: E 2 )^ as the formal difference E\ — 
E 2 . We now define the K-group as the quotient 



K(M) = (V K (M)x V K (M))/~, 



(41) 



which identifies all formal differences that are equiv- 
alent in the sense of Eq 



(|40j). Due to Eq. (|39j), 

every group element in K{Ai) can be represented in 
the form (E,Z n ). However, (E,Z n ) * (E,Z m ) for 
n 7^ m. We define the virtual dimension of (E, F) as 
d v = rk(E) — rk(F), where rk denotes the rank, i.e., 
the fiber dimension of a vector bundle. By restricting 
K(A4) to elements with d v = 0, we obtain the restricted 
K-group K(M) = {g £ K{M)\d v (g) = 0}. K(M) is 
isomorphic to the set of stable equivalence classes 
of Vk{M). Up to now, the construction has been 
independent of the field over which the vector spaces are 
defined. In the following, we will distinguish the real 
and complex K-groups Kr ( M ) , Kq (M ) . 

A crucial notion in K-Theory which is also our main 
physical motivation to study it is that of the stable range. 
The idea is that at sufficiently large fiber dimension n 
no "new" bundles can be discovered by looking at even 
larger fiber dimension. Sufficiently large in terms of the 
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dimension d of A4 means n > nc = d/2 + 1 for the 
complex case and n > = d + 1 for the real case, 
respectively. More formally, every bundle E with n > 
n,K can be expressed as a sum 



E ~ F (B Z r ' 



(42) 



of a bundle F with fiber dimension xik and a trivial bun- 
dle for K = K,C. Since clearly £ ~ F (see Eq. (38 1), 
this means that all stable equivalence classes have rep- 
resentatives in fiber dimension n < n^. Furthermore, a 
situation like our counterexample above where we aug- 
mented two non-isomorphic bundles by the same trivial 
bundle NS 2 to obtain the same trivial bundle cannot oc- 
cur in the stable range. That is to say F as appearing 
in Eq. (42) is uniquely defined up to isomorphisms. The 



stable range hence justifies the approach of if-theory of 
ignoring fiber dimension when defining the stable equiva- 
lence. The key result in the stable range which connects 
K-Theory to our goal of classifying all inequivalent vector 
bundles with sufficiently large but arbitrary fiber dimen- 
sion on equal footing reada^ 



K K (M) = Vect n (M,K) = tt [M,G n (K°°)} 



" (43) 



The complex Bott periodicity Eq. ( 36 1 with period pc = 
2 has a real analogue concerning the homotopy groups of 
O(n) with period pk = 8. This immediately implies in 
the language of K-Theory 



Kk(S' 



K K {S d ), K 



We define 



Kx d (M) = K K {S d M), 



(44) 



(45) 



where S is the reduced suspension^ which for a sphere 
S k indeed satisfies SS k = S k+1 . The stron ger version of 
the Bott periodicity in K-Theory now reads®! 



K 



-d-p K 



K 



(M) 



k-\M), k 



(46) 



which only for A4 = S trivially follows from Eq. 
Using this periodicity, the definition of K^ d in Eq. 
can be formally extended to d € Z. 



The Bott clock 



Class 


Classifying Space 


A 


Co = U(n + m)/ (U(n) xU{m)) 


AIII 


Ci = U(n) 


AI 


R = 0{n + m)/ [0(n) x 0(m)) 


BDI 


Ri = O(n) 


D 


R 2 = U(2n)/U(n) 


Dili 


R 3 = U(2n)/Sp(2n) 


All 


R 4 = Sp(n + m)/ (Sp{n) x Sp{m)) 


CII 


R 5 = Sp(n) 


C 


R 6 = Sp(2n)/U(n) 


CI 


R 7 = U(n)/0(n) 



TABLE III. Table of all classifying spaces C q , R q of complex 
and real K-Theory, respectively. The first column denotes 
the CAZ symmetry class. From top to bottom, the next com- 
plex/real classifying space is the loop space of its predecessor, 
i.e., Oj+i = fiC" 9 (mod 2), R q+1 = flR 9 (mod 8) 



components ir (U(n)) , 7Ti (U(nj) = Tr (£lU(n)) of the 
unitary group and its first loop space, we can classify 
all U(n) vector bundles over S d in the stable range, i.e., 
with n > |. The real analogue of the Bott periodic- 
ity with period pn = 8 leads to analogous statements 
for 0(n) bundles over S d which depend only on the con- 
nected components of 0(n) and its first seven loop spaces 
il l O(n), i = 1, . . . , 7 (see Tab. [ill]). This defines a Bott 
clock with two ticks for the complex case and eight ticks 
for the real case, respectively. Interestingly, these ten 
spaces, for the complex and real case together, are pre- 
cisely the ten Cartan symmetric spaces in which the time 
evolution operators associated with Hamiltonians in the 
ten CAZ classes lie. After this observation, only two 
points are missing until a complete classification of all 
TSM of continuum models can be achieved. The first 
point is a subtlety related to the interdependence of the 
two wave vectors k and — k as shown in Eq. (33), which 



makes the real Bott clock tick counter clockwise. The sec- 
ond point is the inclusion of symmetry constraints into 
the scheme which leads to the clockwise ticking Clifford 
clock (see Eq. ( [50| below). The combination of both im- 
plies that the topological invariant of a continuum model 
of dimension d in the CAZ class q only depends on the 
difference q — d (mod 8) for the eight real classes and on 
q — d (mod 2) for the two complex classes A and AIII, 
respectively. 



From the very basic construction of homotopy groups 
the following identities for the homotopy of a topological 
space X are evident: 

7r d (x) = Tt[s d ,x] = tt [ss d -\x] =tr[s d - 1 ,nx] , 

where Q.X denotes the loop spaced of X, i.e., the space 
of maps from S 1 to X. Iterating this identity gives 
iTd{X) — ir (£l n X) using the complex Bott periodicity 
( 36 1 , we immediately see that counting the connected 



Reality and k-space topology 

For systems which obey anti-unitary symmetries the 
real structure of the Hamiltonian H(x) is most con- 
veniently accounted for in its Majorana representation. 
H(x) can in this representation be expressed in terms of 
a real antisymmetric 2n x 2n-matrix B, 



i 



(47) 
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where c X) ,, i = 1, . . . , 2n are the Majorana operators rep- 
resenting the n fermion species at x. On Fourier trans- 
form, % = J ^ H^S can be written as^ 



U = - J d d fc^'(fc) c _ M c fej , 

where A is skew hermitian and satisfies 
A*(k) = A(-k). 



(48) 



(49) 



Eq. ( 49 ) naturally leads to a real vector bundle structure 



as defined in Ref. [41] for the bundle of eligible ^4-matrices 
over the fc-space (M. d , r), where the involution r (see Ref. 
l4"Tj) is given by k i— > — fc. On one-point compactification, 
this real fc-space becomes a sphere 5 d = (S d ,T) with the 
same involution 17 . Whereas the ordinary sphere S d can 
be viewed as a reduced suspension S of S^ -1 over the 
real axis, S d can be understood as the reduced suspen- 
sion S of S^ 1 over the imaginary axis. This picture is 
algebraically motivated by comparing the involution r to 
the ordinary complex conjugation which, restricted to the 
imaginary axis of the complex plane is of th e sa me form. 



Interestingly, in the lang uage o f definition (45), S plays 
the role of an inverse to < j 17 * 41 *, i.e., 

Kr{M) = K^(SM) = K R (SSM). 

This means that the Bott clock over S is reversed as 
compared to its analogue over S d . 



Real K-theory and the Clifford clock 



base space is the real Brillouin zone T d = (T d ,r), where 
the involution r giving rise to the real structure is again 
given by k h-> — k. For T d the reduced suspension does 
not provide a trivial relation between the K-Thcory of 
different spatial dimension like SS d = S d+1 for the base 
space of continuum models. The general calculation of 
all relevant K-groups over T d has been reported in Ref. 
[T7I Interestingly, the resulting groups always contain the 
respective classification of continuum models in the same 
symmetry class as an additive component. Additionally, 
the topological invariants of weak TSM, i.e., TSM which 
are only present in translation-invariant systems, can be 
inferred. The Clifford clock defined in Eq. (501 is in- 



dependent of the base space and hence still applicable. 
Here, we only review the general result calculated in Ref. 



K R q (f d ) ~ K M q (S a ) 



q( Q d\ 




(51) 



The second term on the right hand side of Eq. (51 ) en- 



tails the notion of so called weak topological insulators 
which are obviously due to TSM in lower dimensions. 
To name the most prominent example, the Z 2 invariant 
characterizing the QSH insulator in d = 2 in the presence 
of TRS, CAZ class All, yields a 3Z 2 topological invari- 
ant characterizing the weak topological insulators with 
the same symmetry in d = 3. 



The main reason for the real construction of Eq. ( 48 ) 



is that the anti-unitary symmetry constraints yielding 
the eight real CAZ classes (all except A and AIII) can 
be distinguished in terms of anti-commutation relations 
of the A-matrix with real Clifford generator d 16 * 17 [ At 
a purely algebraic level, these constraints can be trans- 
formed so as to be expressed only in terms of positive 
Clifford generators^, i.e., generators that square to plus 
identity. We call the restricted K-group of a vector bun- 
dle of ^-matrices over M. that anti-commute with q pos- 
itive Clifford generators K^(A4). Interesting!} * 16 * 17 ! 



ki{M)^k-\M). 



(50) 



Eq. (50) defines a Clifford clock that runs in the oppo- 
site direction as the S d Bott clock. This algebraic phe- 
nomenon explains the full periodic structure of the table 
of TSM of continuum models (see Tab. [TTJ) . The classify- 
ing spaces of A-matrices for systems that anti-commute 
with q Clifford generators are shown in Tab. |III| 



Periodic systems 

The classification of periodic systems is much more 
complicated from a mathematical point of view. Their 



Lattice systems with disorder 

In a continuum model, disorder that is not too short 
ranged so as to keep the fc-space compactification for 
large k valid, can be included into the model system with- 
out changing the classification scheme. However, per- 
turbing a translation-invariant lattice system with disor- 
der also gives its fc-space a discrete lattice structure which 
is not directly amenable to investigation in the framework 
of K- Theory which we only defined over smooth base- 
manifolds. Ref. [T7] shows that a Hamiltonian featuring 
localized states in the energy gap can be transformed 
into a gapped Hamiltonian upon renormalization of pa- 
rameters. The physical consequence of this statement is 
that only a mobility gap is needed for the classification 
of a TSM and no energy gap in the density of states. 
Furthermore, Ref. ^] argues without explicit proof that 
the classification problem of gapped lattice systems with- 
out translation invariance is equivalent to the classifi- 
cation problem of continuum models. This statement 
agrees with the physical intuition that the breaking of 
translation-invariance must remove the additional struc- 
ture of weak TSM as described for periodic systems by 



Eq. (51). 
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VI. CALCULATION OF TOPOLOGICAL 
INVARIANTS OF INDIVIDUAL SYSTEMS 

In Section [Vj we have shown how many different TSM 
can be expected in a given spatial dimension and CAZ 
class. Now, wc outline how insulating systems within 
the same CAZ class and dimension can be assigned a 
topological equivalence class in terms of their adiabatic 
connection defined in Eq. ^ and their adiabatic cur- 
vature defined in Eq. (20 1, respectively. A complete 



case by case study in terms of Dirac Hamiltonian repre- 
sentatives of all universality classes of this problem has 
been reported in Ref. El We outline the general patterns 
relating the classification of neighboring (see T ab. [TT|) 
universality classes following the analysis in RefsP^ZMl, 
Interestingly, all topological invariants can be calculated 
using only complex invariants, namely Chern numbers 
and chiral unitary winding numbers. The anti-unitary 
symmetries are accounted for by the construction of a 



dimensional hierarchy in Section VI B starting from a so 
called parent state in each symmetry class for which the 
complex classification concurs with the real classification. 
In Section [VI C[ we show how the topological invariants 
can be defined for disordered systems with the help of 
twisted boundary conditions. Furthermore, we discuss a 
generalization of the non-interacting topological invari- 
ants to interacting systems in Section [VI D| 



A. Systems without anti-unitary symmetries 

Chern numbers of unitary vector bundles 



Eq. (35) shows that every U{n) bundle E — > M. can 



be represented as a pullback from the universal bun- 
dle £ — > G„(C°°) by some bundle map /. Chern 
classes are de Rham cohomology classes, i.e., topologi- 
cal invariants^! that are defined as the pullback of cer- 
tain cohomology classes of the classifying space G n (C°°). 
The cohomology ring H* (G n (C°°)) consists only of even 
classes and is generated by the single generator Cj € 
H 2 i (G n (C°°)) , j = 1, . . . , n for every even cohomology 
groupP^. The Chern classes Ci of E are defined as the 
pullback Ci — f* 5j from the classifying space by the map 
/ : M — > G n (C°°) associated with the bundle map /. 
Due to the Chern- Weyl theorem^!, Chern classes can be 
expressed in terms of the curvature, i.e., in our case, the 
adiabatic curvature J- defined in Eq. (20) of E. Explic- 
itly, the total Chern class c can be expressed as^ 

< d-l I 1 + =l + c 1 (F)+c 2 (F).... (52) 

The determinant is evaluated in gauge space and prod- 
ucts of T are understood to be wedge products. Cj is 
the monomial of order j in T . Obviously, Cj is a 2j-form 
and can only be non-vanishing for 2j < d, where d is 
the dimension of the base manifold M., i.e., the spatial 



dimension of the physical system. Another characteris- 
tic class which generates all Chern classes is the Chern 
character^! 



ch = Tr 



= l + chi(J")+ch 2 (J') + .. 



(53) 



Due to their importance for later calculations, we ex- 
plicitly spell out the first two Chern characters chi = 
Tr , ch 2 = -giyTr [J 7 A J 7 ]. Importantly, for even 
d — 2p, the integral 



r — 



ch r , 



M 



yields an integer, the so called p-th Chern number^. 
These Chern numbers characterize systems in the uni- 
tary symmetry class A which can only be non-trivial in 
even spatial dimension (see Tab. [TT|. 



Winding numbers of chiral unitary vector bundles 

In Section [VJ we have shown that the classifying space 
for a chiral unitary (AIII) system is given by U(n) and 
that the topological sectors are defined by homotopically 
distinct maps k h-> q(k) € U(n). Now, we discuss how 
to assign a given map q to an equivalence class by cal- 
culating its winding numbe r 44 " 46 ' following Ref. [6[ From 
Tab. [Tl] it is clear that only in odd spatial dimension 
d = 2j — 1 there can be a non-trivial winding number. 
We define 



(-Cj-i)O 



J 2j-1 



(2j - l)!(27ri)J 



^-Tr U'dqf^} , 



(54) 



which has been dubbed winding number density^. Inte- 
grating this density over the odd-dimensional base man- 
ifold M. representing the fc-space of the physical system, 
we get the integral winding number z^j-i 



(55) 



V2j-l = / U>2j-li 
' M 



which is well known to measure the homotopy of the map 
k i y q(k). 



Relation between chiral winding number and Chern Simons 

form 

So far, the relation between the adiabatic connection of 
a chiral system and its topological invariant has not been 
made explicit. Since characteristic classes like Chern 
characters are closed 2j-forms, they can locally be ex- 
pressed as exterior derivatives of a (2j — l)-forms. These 
odd forms are called the Chern Simo ns forms associated 
with the even characteristic clas d 23 * 47 ^. For the j-th Chern 
character chj, which is a 2j form, the associated Chern 
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Simons form Q27-1 rea> 



Q 2j -i(A,? t ) = 



U - 1)! 



dt STr 



A, J-j 



(56) 



where Ft = t T + (t 2 — t)A A A is the curvature of the 
interpolation tA between the zero connection and A and 
STr denotes the symmetrized trace. Explicitly, we have 
Q 1 = ^Tr [A] , Q 3 - - g^Tr [AdA + f „4 3 ] . 

It is straightforward to show^, that in a suitable gauge, 
the Berry connection of a chiral bundle yields A B = 
\qdq* , where q € U(n) is again the chiral map char- 
acterizing the system. This is not a pure gauge due to 
the factor \ which entails that the associated curvature 



T B does not vanish. Plugging A B and JF B into Eq. ( 56 ) 
immediately yields^ 



Q2j-i(A B ,T B ) 



B\ q 



(57) 



Eq. (57) directly relates the winding number density to 
the Chern Simons form. We define the Chern Simons 
invariant of an odd dimensional system as 



CS2J-1 



M 



Q 2j -i (mod 1), 



where (mod 1) accounts for the fact that J M Q.2j~i has 
an integer gauge dependence due to i^2j-\ (U(n)) = Z for 
n > j. Looking back at Eq. (55), we immediately get 



j/ 2 j_i(mod 2) = 2CS 2 j_i(mod 2). 

We note that the (mod 2) can be dropped if we fix the 
gauge as described above to A B = ^qdq^ 1 . This estab- 
lishes the desired relation between the winding number 
of a chiral unitary system and its adiabatic curvature. 



B. Dimensional hierarchy and real symmetry 
classes 



Until now, we have only discussed how to calculate 
topological invariants of systems in the complex sym- 
metry classes A and AIII. Interestingly, for some real 
universality classes, the classification in the presence of 
anti-unitary symmetries concurs with the unitary classi- 
fication (see Tab. IV). The first known example of this 
type is in the symplectic class All in d = 4 which is char- 
acterized by the se cond C hern number of the correspond- 
ing complex bundle^^. Another example of this kind is 
the p -Hp superconductor in d — 2 and symmetry class D 
which is characterized by its first Chern number, i.e., in 
the same way as the QAH effect in class A. In odd dimen- 
sions similar examples exist for real chiral classes, e.g., for 
Dili in d = 3, where the winding number is calculated 
using Eq. ( 54 ) in the same way as for the chiral unitary 



v = 




TRS 



FIG. 1. Illustration of the WZW dimensional extension. The 
circle at v = represents the physical system. The poles 
at v = ±7r represent the trivial reference system without k- 
dependence. The two interpolations are conjugated by an 
anti-unitary symmetry, here exemplary denoted by TRS. 



class AIII in the same dimension. All the topological in- 
variants just mentioned are integer invariants. In some 
universality classes, these integers can only assume even 
values (see Tab. IV). For physically relevant dimensions, 
i.e., d = 1,2,3, these exceptions are CII in d = 1, C in 
d = 2, and CI in d — 3. All other states where the com- 
plex and the real classification concur, can be viewed as 
parent states of a dimensional hierarchy within the same 
symmetry class from which all Z 2 invariants appearing 
in Tab. |IV| can be obtained by dimensional reduction. 
This approach was pioneered in the seminal work by Qi, 
Hughes, and ZhangSZl. 



The general idea is more intuitive if we consider 
the parent state as a W ess-Zumino-Wittcn (WZW) 
dimensional extensio n 1 50 * 51 ! of the lower dimensional 
descendants instead of thinking of a dimensional re- 
duction from the d-dimensional parent state to its 
descendants. This works as follows: We fix a localized 
(d — l)-dimensional insulator without any hopping 
that satisfies the required anti-unitary symmetries 
as a trivial reference state. This reference state is 
described by the independent Bloch Hamiltonian ho. 
The (d — l)-dimensional physical system of interest is 
characterized by the Bloch Hamiltonian h(k). Then, 
we interpolate by varying the parameter v between 
the (d — l)-dimensional physical system (v = 0) and 
the trivial state (v — 7r) without closing the insulating 
gap. However, the intermediate (d — l)-dimensional 
system at fixed v 7^ 0, 7r might well break the required 
anti-unitary symmetries. The crucial step is now to 
do the interpolation for v G [0, tt] and v £ [—it, 0] in a 
symmetry conjugated way (see Fig. [T]). That is to say, 
we require our (d — l)-dimensional system of interest 
and the resulting <i-dimensional extended system to be 
in the same CAZ class. This d-dimcnsional system is 
characterized by the Bloch Hamiltonian h(k,v). The 
v € [— 7r,0] and the v £ [0,7r] half of the extended 
fc-space then are not independent of each other but give 
equal contributions to the integer top olog ical invariant 
of the d-dimensional extended systenPEU. One might 
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Class 
A 

AIII 



AI 
BDI 

D 

Dili 

All 
CII 

c 

CI 



d=ld=2d=3d=4d=5d=6d=7d=8 







z 









z 





2Z 





z 2 
z 2 



2Z 




2Z 






^2 

z 2 



2Z 





2Z 






z 2 





z 2 



2Z 






z 2 
z 2 



1L_ 

z 2 
z 2 



2Z 






TABLE IV. Table of all groups of topological equivalence classes. The first column denotes the symmetry class, divided into 
two complex classes without any anti-unitary symmetry (top) and eight real classes with at least one anti-unitary symmetry 
(bottom). Chiral classes are denoted by bold letters. The parent states of dimensional hierarchies are boxed. For all non-chiral 
boxed states, the classification concurs with that of class A in the same dimension. For all chiral boxed states, the classification 
concurs with that of class AIII in the same dimension. 2Z indicates that the topological integer can only assume even values 
in some cases. Such states are never parent states. 



now ask to which extend the resulting integer invariant 
of the extended system depends on our choice of the 
interpolation h(k,v) between h(k) — h(k,v — 0) and 
ho = h{k,v — ±7r). To answer this question, one 
considers two interpolations h(k,v),h(k,v). It is then 
elementary to show^ that the difference between the 
integer invariants of these two d-dimensional systems 
is an even integer. This implies that a Z2 information, 
namely the parity of the integer invariant associated 
with the extended system, is well defined only in terms 
of the physical system with spatial dimension [d — 1). 

A similar procedure can be repeated a second time to 
obtain a Z2 invariant for a (d — 2) dimensional second 
descendant!^. From the procedure just sketched, it is 
obvious why the exceptional phases which are character- 
ized by an even integer are not parent states of such a 
dimensional hierarchy. The generic const ructions for all 
possible classes can be found in RefsP 23 . With that, we 
are provided with a general and fairly explicit recipe for 
the practical calculation of the topological invariants for 
all possible CAZ classes in all spatial dimensions. 



C. Bulk invariants of disordered systems and 
twisted boundary conditions 

Our practical calculation of topological invariants 
so far has been focused on periodic systems with a 
BZ, i.e., A4 = T d and continuum models where the 
/c-space can be compactified to a sphere, i.e., A4 = S d . 
As already pointed out in Section |Vj the situation 
is more complicated for disordered lattice models. 
Seminal progress along these lines was reported for the 
quantum Hall state by Niu et al. in 1985^1. These 
authors use twisted boundary conditions (TBC) to 
define the quantized Hall conductivity a xy for a 2D 



system as a topological invariant only requiring a 
bulk mobility gap. We briefly review their analysis 
and propose the framework of TBC as a general recipe 
to calculate topological invariants for disordered systems. 

The Hall conductivity a xy resulting from a linear re- 
sponse calculation at zero temperature yields 

_-2 \- (Ol^lnXnl^lO) 

where |0) is the many body ground state, A is the area 
of the system, and Hi = §jj~. In the presence of a mag- 
netic field, translation-invariance is defined in terms of 
the magnetic translation operator tJ^Q which concurs 
with the ordinary translation operator T(a) = e lak in 
the absence of a magnetic field. TBC now simply mean 
that a (magnetic) translation by the system length Lj in 
j-direction gives an additional phase factor e l( ^ j . <pj is 
called the twisting angle in j-direction. Gauging away 
this additional phase to obtain a wave function with pe- 
riodic boundaries amounts to a gauge transformation of 
the Hamiltonian which shifts the momentum operator 
like 

— idj —> —idj + <j)j (58) 

Using Hi = Hfc = [dfafH] the Hall conductivity can 
be expressed as the sensitivity of the ground state wave 
function to TBC. 

2 

a xy = -Im(9^O|5 0y O), 

or, defining 9 = L x <f> x , <p = L y ^> y 

a xy = 2Im(d e 0|cV0) = uF 0v . (59) 

The main merit of Ref. [53] is to show that this expression 
actually does not depend on the value of (9, <p) as long 
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as the single particle Green's function of the system is 
exponentially decaying in real space. This condition is 
met if the Fermi energy lies in a mobility gap. Hence, a 
trivial integration can be introduced as follows: 



extension (see Section VI B I , a hybrid approach between 



i 

4tt2 



2tt 



2tt 



f iT 

dfJ 7 e v = G — = G Ci, 
J T 2 Ztt 

(60) 



2 -I 

where Go = %- = ^ ancl the integer C\ is by definition 
the first Chern number of the ground state line bundle 
over the torus of twisting angles. This construction 
makes the topological quantization of the Hall conduc- 
tivity manifest. 



Eq. ( 58 ) shows the close relation between momentum 



and twisting angles. One is thus tempted to just replace 
the BZ of each periodic system by the torus of twisting 
angles for the corresponding disordered system which is 
topologically equivalent to a fictitious perio dic sy stem 
with the physical system as single lattice sit d 53 * 54 *. We 
will proceed along these lines below but would like to 
comment briefly on the special role played by the quan- 
tum Hall phase first. Eq. (59) represents the physical 



observable a xy in terms of the twisting angles. Niu et al. 
argued rigorously^! that a xy of a bulk insulating system 
can actually not depend on the value of these twisting 
angles which allows them to express a xy as a manifestly 
quantized topological invariant in Eq. ( |60| . For a 
generic TSM, the topological invariant of the clean sys- 
tem does in general not represent a physical observable. 
Furthermore, the integration over the twisting angles 
will not be trivial, i.e., the function to be integrated will 
actually depend on the twisting angles. Employing the 
picture of a periodic system with the physical system 
as a single site is problematic inasmuch as the bulk 
boundary correspondence at the "boundary" of a single 
site is hard to define mathematically rigorously. In the 
quantum Hall regime for example, it is well known that 
in disordered systems a complicated landscape of local- 
ized states and current carrying regions produces the 
unchanged topologically quantized Hall conductivity®. 
However, the edge states of the disordered quantum Hall 
state are in general not strictly localized at the boundary. 

Replacing the BZ of a translation-invariant system by 
the torus of twisting angles in the disordered case yields 
a well defined topological invariant which adiabatically 
connects to the topological invariant of the clean system 
where the relation 



(61) 



follows from Eq. ( 58 1 . This is because from a purely 



mathematical perspective it cannot matter which torus 
we consider as the base space of our system. In this 
sense, the framework of TBC is as good as it gets con- 
cerning the definition of topological invariants for disor- 
dered systems. When calculating the topological invari- 
ant of a symmetry protected TSM through dimensional 



TBC in the physical dimensions and periodic boundaries 
in the extra dimension can be employed to define a topo- 
logical invariant for the disordered system which can be 
calculated more efficientlySSl. The fact that in some sym- 
metry protected topological phases the topological invari- 
ants are not directly representing physical observables is 
a not a problem of the approach of TBC but is a remark- 
able difference between these TSM and the quantum Hall 
state at a more fundamental level. 



D. Taking into account interactions 



Up to now, our discussion has only been concerned 
with non-interacting systems. As a matter of fact, 
the entire classification scheme discussed in Sec. |V] 
massively relies on the prerequisite that the Hamiltonian 
is a quadratic form in the field operators. The violation 
of this classification scheme for systems with two particle 
interactions has been explicitly demonstrated in Ref. [56l 

As we are not able to give a general classification of 
TSM for interacting systems, we search for an adiabatic 
continuation of the non-interacting topological invariants 
to interacting systems. This procedure does from its 
outset impose certain adiabaticity constraints on the 
interactions that can be taken into account. The topo- 
logical invariants for non-interacting systems are defined 
in terms of the projection P on the occupied single 
particle states defining the ground state of the system. 
The main assumption is thus that the gapped ground 
state of the non-interacting system is adiabatically 
connected to the gapped ground state of the interacting 
system. A counter-example o f this phenomenology are 
fractional quantum Hall state^^H, where a gap due to 
non-adiabatic interactions emerges in a system which is 
gapless without interactions. However, it is clear that 
the phase space for low energy interactions will be much 
larger in a gapless than in a gapped non-interacting 
system. We thus generically expect the classification 
scheme at hand to be robust against moderate inter- 
actions. However, beyond mean field interactions, the 
Hamiltonian cannot be expressed as an effective single 
particle operator. Hence, we need to find a formulation 
of the topological invariants that adiabatically connects 
to the non-interacting language and is well defined 
for general gapped interacting systems. The key to 
achieving this goal is to look at the single particle 
Green's function G instead of the Hamiltonian. This 
approach has been pioneered in the field of TSM by Qi, 
Hughes, and ZhangS3 who formulated a topological field 
theory for TSM in the CAZ class AIL 
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Chern numbers and Green's function winding numbers 

The role model for this construction is again the Hall 
conductivity a xy of a gapped 2D system. In Ref. |44j 
a xy has been expressed in terms of G by perturbative 
expansion of the effective action of a gauge field A that 
is coupled to the gapped fermionic system. The lead- 
ing contribution stemming from a vacuum polarization 
diagram yields the Chern Simons action 



d 2 xdt A„d v A a = -2L I A ,IA. 



The prefactor a xy in units of the quantum of conductance 
assumes the form 



'■iy 



2 l — 2 J d 2 fcdcTr [(GdG- 1 ) 3 } , (62) 



where d now denotes the exterior derivative in combined 
frequency-momentum space and G is the time ordered 
Green's function, or, equivalently as far as the calcula- 
tion of topological invariants is concerned, the continu- 
ous imaginary frequency Green's function as used in zero 
temperature perturbation theory. An analogous expres- 
sion has been derived by Volovik using a semi-clas sica l 
gradient expansion^. The similarity between Eq. (54) 



and the integrand of Eq. (62) is striking. Obviously, Eq 



( 62 ) represents a xy as a winding number in 3D frequency- 
momentum space. If this construction makes sense, we 



should by integration of Eq. ( 62 ) over uj recover the rep- 



resentation of a xy as the first Chern number in the 2D 
BZ for the special case of the non-interacting Green's 
function Go(uj,k) — (ius — h(k))~ 1 . This straightforward 
calculation relies on the residue theorem and has been 
explicitly presented in Ref. [37| Its result can be readily 
generalized to higher even spatial dimensions and higher 
Chern numbers, respectively. In Ref. a perturba- 
tive expansion similar to Ref. 2J] has been presented 
for fermions coupled to a gauge field in arbitrary even 
spatial dimension 2n. The resulting analogue of the Hall 
conductivity, i.e., the p refa ctor of the Chern Si mons form 
in (2n + 1)D (see Eq. {56])) can be expressed a d 10 ' 37 ' 46 ' 



2r+l 



[G] = N (r) 



Tr 



BZ> 



GdG- 



(63) 



Af(r) = 



-r-l 



(2r + l)!(27ri)' r + 1 ' 



Performing again the integration over the frequency ana- 
lytically for the noninteracting Green's function Go yields 



A^r+l [Gq] — C r 



(64) 



Eq. (64) makes manifest that A^ r +i [G], which 



can be formulated for an interacting system, repro- 
duces the non-interacting classification for the free 
Green's function Go of the non-interacting system. 
The topological invariance of N2 r +i [G] is clear by 



analogy with Eq. (55): Whereas the winding num- 



ber V2j-i measures the homotopy of the chiral map 
k n> q(k) 6 U (n) which, properly normalized, yields 
an integer due to ^2j-i {U(n)) — Z, n > j, Eq. (63) 
measures the homotopy of G G GL(n + m, C) in the 
(2r + 1)D frequency-momentum space which is also 
integer due to 7r 2r +i (GL(n + m, C)) = Z, n + m > r. 

The dimensional hierarchy for symmetry protected de- 
scendants of a parent state which is characterized by a 
Chern number (see Section VI B ) can be constructed in 
a completely analogous way for the interacting general- 
ization N2r+i of the Chern number Cj^. The resulting 
topological invariants for the descendant states have been 
coined topological order parameters in Ref. 1511 Disor- 
der can again be accounted for by imposing TBC and 
replacing the fc-space of the system by the torus of twist- 
ing angles (see Section VI C I . 



Interacting chiral systems 

The integer invariant of chiral unitary systems (class 
AIII) in odd spatial dimension 2r — 1 is not a Chern 



number but a winding number (see Section VI A ) . For all 
these systems and dimensional hierarchies with a chiral 
parent state, i.e., all chiral TSM (see Tab. IV), a similar 
interacting extension of the definition of the invariants in 
terms of G(iui, k) has been reported in Ref. [SHI 



hr [G] = n(r) 



Tr 



BZ> 



Q(dQ) 



2r 



(65) 



where n(r) is a normalization constant, and Q(iu), k) — 
G^ 1 (ioj)UciiG(iu!,k) 1 with the unitary representation 
matrix Uqh of the chiral symmetry operation. In the 
non- interacting limit, l2 r reduces to f2r—i as defined in 
Eq. ((55pl. 



Fluctuation driven topological transitions 

Thus far, we have shown that for a non-interacting 
system the integration over uj reproduces the band struc- 
ture classification scheme formulated in terms of the adi- 
abatic curvature. However, the additional frequency de- 
pendence of the single particle Green's function can cause 
phenomena without non-interacting counterpart. To see 
this, we represent the single particle Green's function of 
an interacting system as 

G(oj, k) = {ioj - h(k) - E(w, A:)) -1 , 

where S is the self-energy of the interacting system. In 
Ref. 59, it has been pointed out, that the value of 
N2r+i [G] cannot only change due to gap closings in the 
energy spectrum as in the case of the Chern number C r . 
This is due to the possibility of poles in the w-dependence 
of the self-energy which give rise to zeros of the Green's 
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function, whereas gap closings correspond to poles of G. 
From the analytical form of N\ it is immediately clear 
that both poles and zeros of G can change the value of N\ . 
More generally, the G <-> G -1 symmetry of N2 r +i makes 
clear that poles of G can be seen as zeros of G -1 and 
vice versa on an equal footing. In Ref. [6Ql it has been 
demonstrated that the w-dependence of S can change a 
non-trivial winding number into a trivial one. The emer- 
gence of a topologically nontrivial phase due to dynami- 
cal fluctuations which has no non-interacting counterpart 
has been presented in Ref. I6T1 



Chern numbers of effective single particle Hamiltonians 

Due to its additional w-integration, the practical cal- 
culation of iVar+i can be numerically very challenging 
once the single particle Green's function of the interact- 
ing system has been calculated. A major breakthrough 
along these lines has been reported in Ref. These 
authors show using the spectral representation of the 
Green's function that one can always get rid of the in- 
dependence of G. We only review the physical results 
of this analysis. The both accessible and explicit proof 
can be found in Ref. The physical conclusion is as 
elegant as simple: Instead of calculating N2 r +i we can 
just calculate the Chern number C r associated with the 
fictitious Hamiltonian 



h(k) = -G _1 (0,A), 



(66) 



the occupied states of which are just its eigenstates with 
negative eigenvalues which have been dubbed i?-zeros^l 
since they are positive energy eigenstates of G _1 (0, k). 
Obviously, — G ( ^" 1 (0, k) — h(k) for the non-interacting 
Green's function. Hence, h(k), which has recently been 
coined topological HamiltoniarP^ adiabatically connects 
to h{k) in the non-interacting limit. Note that the pos- 
sibility of eliminating the w-dependence is not in con- 
tradiction to the relevance of this ^-dependence for the 
topology of the interacting system. All it shows is that 
the relevant changes due to a different pole structure of 
G as a function of u can be inferred from its value at 
u> = 0. 



Topological Hamiltonian for chiral interacting systems 

In principle, the construction of the topological Hamil- 
tonian h(k) can be readily generalized to chiral interact- 
ing systems as has been mentioned in Ref. H>U To see 
this, we note that the crucial argument for the construc- 
tion of the topological Hamiltonian brought forward in 
Ref. (62] is the following: The continuous interpolation 

G(iuj,k,X) = (l-\)G(iuj,k) + \[iuj + G- 1 (0,k)]~ 1 

does not contain any singularities or gap closings. Thus, 
as long as the calculation of a topological invariant 



in terms of G(iu>, k) is concerned 
G(iu),k) = G(iw,k,X = 1) 



viously, G(iu), k) — 



h{k) 



-l 



we can also use 
G-^fc)] -1 . 



Ob 

is the Green's func- 



tion of a fictitious non-interacting system which is gov- 
erned by the topological Hamiltonian h(k). The mere 
existence of the topological invariant for chiral systems 
in terms of G(iui, k) as presented in Ref. [SH (see also Eq. 
(65 1) hence suffices to argue that one can equally well 



investigate the topology of h(k) and its symmetry pro- 
tected descendants (see Section |VIB ) instead of directly 
evaluating Eq. ( [65] ) . Since the single particle Green's 
function inherits the fundamental symmetries from the 
HamiltoniarF^I, h(k) will also obey these symmetries. In 
particular, for an interacting system with chiral symme- 
try, the topological Hamiltonian can be brought into the 
fiat band off-diagonal form 



h(k) 



q(k) 
fl~t(fc) , 



(67) 



where q(k) E U(n) for the topologically equivalent 
flat-band system. This construction generically defines 
a topological invariant for the chiral interacting system 
which adiabatically concurs with the non-interacting 
system: The winding number v associated with the 
fictitious Hamiltonian h{k). A similar construction for a 
chiral ID system has been presented in Ref. [551 



Discussion of the topological Hamiltonian and practical 
consequences 

As already mentioned, the above construction cannot 
be valid for arbitrary interacting systems. In particular 
in ID, the breakdown of the Z classification in the 
presence of general interactions has been investigated 
in Ref. 1561 However, this problem does not pertain 
to the concept of the topological Hamiltonian itself 
but rather reflects the limited validity of the adiabatic 
continuation of the non-interacting invariants in terms 
of the single particle Green's function, i.e., the limited 
validity of Eq. (63) and Eq. (65). In the validity 



regime of these equations, one can cquivalcntly use the 
topological Hamiltonian h(k) = — G _1 (0, k) to classify 
an interacting system in any symmetry class. This is of 
enormous practical usefulness for at least two reasons. 
First, we get rid of the w-integration appearing in Eq. 
(63) and Eq. (65) which is cumbersome to evaluate. 



Second, the method of dimensional extension, though 
generally valid, is not always the most convenient way 
to calculate the topological invariant of a symmetry 
protected descendant state. Provided with the formal 
equivalence between the non-interacting classification 
problem of the topological Hamiltonian and the Green's 
function topology, we can directly apply all simplified 
schemes that have been introduced to directly calculate 
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non-interacting invarian ts of symmetry protected states 
(see, e.g., Refs! 31 ' 53 1 66 ^ to the topological Hamilto- 
nian. The framework of dimensional extension and Eq. 
(63) or Eq. (65) for the parent state are, with the benefit 



of hindsight, only needed to justify the validity of the 
topological Hamiltonian. 

Before closing the section, we would like to discuss the 
role of the bulk boundary correspondence, in the pres- 
ence of interactions. In general, interactions can spon- 
taneously break the protecting symmetry of a symmetry 
protected TSM locally at the boundary thus gapping out 
the characteristic metallic surface states. Importantly, 
this spontaneous symmetry breaking at the gapless sur- 
face will typically happen at a lower critical interaction 
strength than in the gapped bulk. This is because the 
gapless surface modes offer more phase space for interac- 
tions. A generally valid bulk boundary correspondence is 
hence absent in the interacting case. Within the validity 



regime of Eq. (63) and Eq. (65) for chiral TSM, respec 



tively, an interacting analogue of the bulk boundary cor- 
respondence has been reported in Refsf^EH. The main 
difference to the non- interacting case is that boundary 
zero-modes, which represent poles of the Green's function 
can be canceled by zeros of the Green's function as far 
as the calculation of topological invariants is concerned. 
Note that the Green's function of a non-interacting sys- 
tem does not have zeros. 



broader context, we would like to point out that there are 
also phenomena of topological origin which emerge only 
due to the presence of interaction s. Th e historically first 
phenomenon is the £ FQH effeciP^l which cannot be 
adiabatically connected to an insulating non-interacting 
state. The non-interacting state is in this case a partially 
filled Landau level which provides an enormous density 
of states at the Fermi energy. In a system with periodic 
boundaries, the £ FQH has a characteristic ^-fold ground 
state degeneracy. Interestingly, this simplest FQH state 
can still be analyzed in the framework of TBCP^l For 
more general FQH systems the concept of topological or- 
der has been introduced by WerPH. A crucial notion in 
this framework is the quantum dimension of the topolog- 
ically ordered system which can be viewed as the ground 
state degeneracy of the system on a torus, i.e., with pe- 
riodic boundary conditions. Most TSM discussed in this 
article have quantum dimension one, like a trivial insu- 
lator. From the vantage point of topological order, these 
states are thus trivial. Recently, the concepts of TSM 
and FQH physics have been combined to the definition 
of the fractional Chern insulator^HZH anc j the fractional 
topological insulatoi^HZSl^ a translation-invariant realiza- 
tion of the FQH and its symmetry protected analogues, 
respectively. 



VII. LIMITATIONS OF THE FRAMEWORK OF 

TSM 

Finally, we would like to point out some limitations 
of the concept of TSM. The two main aspects that one 
could see critical in the field of TSM are outlined in the 
following. 

First, whereas the topologically quantized Hall 
conductivity in the integer quantum Hall state, the 
historical role model of all TSM, is a physical observable, 
the topological invariants of symmetry protected TSM 
like the QSH state are not directly physically observable 
without additional unitary symmetries. The quantum 
Hall effect can be understood in terms of the spectral 
flow associated with the threading of a flux tube^l. 
Along similar lines, the QSH effect can be understood 
in terms of a spin charge separation associated with 
the threading of a spin flru>p2. However, this spin flux, 
as opposed to an ordinary magnetic flux tube, has no 
known experimental realization. 

Second, the entire construction and classification of 
TSM is based on single particle Hamiltonians. In Sec- 
tion VI D we discussed how adiabatic interactions can be 



taken into account and argued that interactions of mod- 
erate strength are not likely to destroy the phenomenol- 
ogy of TSM. In order to position the field of TSM in a 



VIII. OUTLOOK 



From a conceptual point of view, the entire zoo of 
topological states of matter can be seen as conclusively 
understood in the framework outlined in this review. 
However, there are at least two general routes to be 
considerably further explored by future research. First, 
the precise experimental implications of many topolog- 
ical states of matter have not been fully unraveled yet. 
Whereas in the quantum Hall state, the topological in- 
variant directly represents a physical observable, namely 
the Hall conductivity of the sample, the observability of 
the topological invariants of several topological states 
of matter is unknown or still under active debate. 
This issue is from our point of view closely related to 
the rather limited number of promising proposals for 
concrete technological applications based on these novel 
states of matter. Obviously, successful research in this 
direction will be of decisive importance for the long 
term future of the entire field of topological states of 
matter. Second, the influence of interactions and open 
quantum system effects on topological states of matter 
is by no means conclusively understood, let alone an 
exhaustive topological classification of interacting or 
dissipative systems. As a first step along these lines, 
a purely dissipation driven topological state has been 
reported in Ref. [BUI 
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